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Abstract 

We derive coupled nonlinear Schrodinger equation (CNLSE) for arbi- 
trary polarized light propagation in a single-mode fiber. We introduce 
a basis of transverse eigen modes with the appropriate projecting hence 
solutions depend on the waveguide geometry. Considering a weak non- 
linearity which is connected with Kerr effect, we give explicit expressions 
for nonlinear constants via integrals of Bessel functions. We compare 
numerical results for the nonlinear constant extracted from experimental 
observations of a soliton for the nonlinear Schrodinger equation (NLSE) 
(single-mode one). The method of projecting we use allows a direct gen- 
eralization to multi-mode fiber case. 



1 Introduction 

There are lot of publications devoted to the propagation and interaction of po- 
larized electromagnetic wave pulses in optical fibers (see the book ). Almost 
most of them exploit the results of [2], that is going up to [3]. It is claimed that 
one can consider the fiber as made of isotropic material and the birefringence is 
originated from the third order nonlinearity (Kerr effect 4 ) . The result of the 
derivation is achieved by means of averaging across the fiber section and gives 
the following evolution along z axis. 

i X+-ik'+X++ , ^X+ + (<y\X+\ 2 + T 1 \X-\ 2 )X+ + ... = 0, (la) 

iX- -ik'^Xr + yA t 7 + ( 7 |A-| 2 + V \X+\ 2 )X- + ... = 0, (lb) 

where X + , X~ are the envelopes components of electric fields (polarizations), 
k" is the dispersion constants respectively, 7 correspond to SPM (self phase 
modulation) and 77 corresponds to XPM (cross phase modulation). The compu- 
tation of the 7/77 relation (ratio) for Kerr medium which is generally elliptically 
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bircfringent, depends on birefringent ellipse axis choice [3 El E] and have 
value 2/3 (linear case), 2 (circular case), generally (for example 7/77 = 2) we 
can say that the XPM is twice as effective as SPM . V^ N = 1/ k'^ is the nonlinear 
group velocity of polarization components. Here we accept that the origin of a 
birefringence comes from nonlinear effects, but it can be also descent from ran- 
dom defects of fiber or special structure of waveguide (polarization maintaining 
fibers). 

The averaging procedure looks reasonable from a physical scope but in many 
cases leads to significant deviations from experiments The transition from 
a three-dimensional to one-dimensional picture by the averaging is quite impos- 
sible in the case of multi-mode field: it leads to the only equation while the 
modes should be described by independent variables. 

In this paper we base on a projecting procedure to the mode subspaces in 
a functional space of a multi-mode field [H]- In the nonlinear theory it leads 
to the important difference between results for nonlinear constants obtained by 
the projecting and averaging procedures already in the case of one-mode fiber. 

The general plan of the paper is following. In the section [21 we briefly 
overview (fix) the notation that are chosen maximally close to standard books 
4 . In the next section |31 we show how to build the general representation of 
the overall electromagnetic field as the mode superposition. We also define the 
transverse orthogonal eigenmodes together with the projecting procedure by 
the appropriate scalar product. The section 01 contains numerical results for a 
single mode waveguide and the evaluation of nonlinear coefficients for the NLS 
equation. 



2 Basic equations 

We describe polarization modes interaction in the cylindrical optical fibers. We 
start from the Maxwell electromagnetic field equations 

VB = 0, (2a) 
V D = 0, (2b) 

VxE = --, (2c) 
<9D 

VxH = — , (2d) 

in the system in the cylindrical polar coordinate (r, tp, z) 
and materials equations 

H = — B, (3a) 
Mo 

D = e E + P. (3b) 

When one study boundary conditions the polarization vector P is considered as 
a linear function of E, we take the simplest form for the isotropic medium 

P = ^OXimearE (4) 
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and a wave equation for the electric field in fiber is 



AE - = 0. (5) 

Boundary conditions for our waveguide are 

D rl -D r2 = 0, (6a) 

B rl -B r2 = 0, (6b) 

n x (Ea - Ex) = 0, (6c) 

nx(H a -Hi) = 0. (6d) 

Conditions l(6"c|) for electric field yield 

e 1 E r (r a+ ) = e 2 E r (r ..), (7a) 

E <p (r 0+ ) = E v (r Q -), (7b) 

E z (rQ +1 ip,z) = E z (ro-,ip,z) (tq - waveguide radius). (7c) 

Wave number k must be the same inside and outside a waveguide. To 
perform boundary conditions, we defined two parameters a and [3 

a 2 = a> 2 eoA i o£i — k 2 , r < ro, (8a) 

f3 2 = k 2 - oj 2 e ^e 2 , r > r , (8b) 

where to is the frequency of a light wave. 



Now if we use solution inside and outside then waveguide for linear polariza- 
tion and all of boundary conditions we get equation (known as Hondros-Debye 
equation) for the eigenvalues ai n (see equations |HJ) . From this equation (which 
is well known in linear theory of waveguides U) we can numerically evaluate 
eigenvalues. The linearized ME defined the basis of eigenfunctions that are 

Ji(a ln r)e llv (9) 

where i = 0,±l,±2,... and n — 1,2,... where n is numbering following eigen- 
values (following solutions for fixed I) and ai n is connected with eigenvalues 
kin ■ 

In general the polarization vector should be written as 

P = e (^X (1) E + x (2) :EE + x (3) :EEE + • • • ^ , (10) 

where \^ is linear dielectric susceptibility and corresponds to the refraction of 
light. In a case of the second order dielectric susceptibility x^ 2 \ we could omit 
it because it is equal zero in materials construct with symmetrical molecule. 
From higher order dielectric susceptibility we save only third order susceptibility 
because rest of orders are negligible. The third order susceptibility is responsible 
for nonlinear refraction of light, self phase modulation (SPM) and cross phase 
modulation (XPM). 
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If we take into consideration that for impulses longer then O.lps one can 
treat a response of a medium as instantaneous and we can write 



PnlW =£oX (3) (M,^E(t)E(t)E(i). (11) 

The third order dielectric susceptibility \^ f° r isotropic media is discussed in 
the papers [31 HH- Basing on it we write 

"X,ijkl — Xxxxx$ij$kl ~\~ Xxyxy^ik^jl H~ "X_xyyx$il$kj •, (^-^) 
\xxxx = Xyyyy = Xzzzz = \xxyy H~ X^yzy "I - \xyyx-, 

(13) 

all components of electric field are in standard form 

Ei = ^A ie iut + c.c, (15) 

inserting this relation into equation 1|11|) , we get nonlinear polarization as (non- 
resonant terms are removed) 

Pi = Ixxxx^oYl ( 2A i\Aj\ 2 + AjAi) + c.c, (16) 

j 

where i,j — x,y,z 

For example, for z component we have 
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Pz T^Xxxxx^-O 



|V + |(|^-| 9 + I^| 9 ) 



A z + l -A z (Al + ADje^ + cc. 

(17) 

We have same equation as in paper |3] [2] . Next we put it into the Maxwell 
equations (We introduce nonlinearity into the Maxwell equations in the form of 
the Kerr effect pQ, with assumption of small nonlinearity.) 
Let us rewrite the wave equation system as 

d 2 

where □ is defined by 

d 2 

□ = A'o£o£^2 - A - ( 19 ) 
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3 General solution, main results 

We write solution for electromagnetic field with amplitude depend on time and 
propagation coordinate in form jS] : 



(20a) 



E r (r, tp,z,t) 



p,l,n ™ " 



^ n (^ t )~ J l( a nl r ) + Cf n (z,t)kd r Ji(a n ir) 



(20b) 



E<p{r,<p,z,t) 



2 V 



iZ A; 



T>f n (z,t)ujd r Ji{a n ir) - £f (z,t) — Ji(a ni r) 



e ilV e i("t-kz) + cc 



(20c) 



B,(r, = \ [^{z^Jiiamry^e^-^ +c.c. 



(20d) 



B r (r,ip,z,t) 
2^ 15 



^ ^ ilcjfx o e e ^ ^ _ ^ ^ t)kd r Ji (a n ir) 



Jltp e i(ut-kz) , ( 



(20e) 



B v (r,ip,z,t) 



iZfc 

7^ (z, t)uj^ o e £d r Ji (a n ir) + Sf n (z, t) — Ji (a„ir) 



e aip e i(u,t-kz) + cc 



(20f) 



Here p numbering two orthogonal polarization and have values "+" and "— ". 
Coefficients with tilde includes all constants to simplify notation. 
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Inserting these solutions into the Maxwell equations yields 





= o, 


(21a) 


Cf n al + d z Al 


= o, 


(21b) 


-cii 2 -sf n ii 


= o, 


(21c) 




= o, 


(21d) 


-d z Vl + d t Ql 


-o, 


(21e) 


- dzC f n -Ai-d t vr n 


= o, 


(21f) 


-d z B? n - d t Sl 


= o, 


(21g) 


-Vf n al + d t T[ n 


= o, 


(21h) 




= o, 


(21i) 


-zi + Ml 


= o, 


(21j) 




= o, 


(21k) 


Hlal+d z T? n 


= o, 


(211) 


- Sfjl 


= o, 


(21m) 


F? n il + d z Sf n + fi e ed t Bf n 


= 0, 


(21n) 


d z Vf n + n e ed t C P n 


-o, 


(21o) 


d z Ql + d t Ef n 


= 0, 


(21p) 


- d *n p ln - T p n - n e ed t V p ln 


= 0, 


(21q) 


afnK t - ^ed t A p ln 


= 0, 


(21r) 


-W 2 - s? n u 


= 0, 


(21s) 




= 0. 


(21t) 



Using (|21b|l . Ij21r|) and (f 2 1 ff) we can verify that the amplitudes Af satisfy 
equation (similarly we can get equation for J-f n ) 

d zz Al-noeoed tt Al = alAl (22) 

It can be proved that the Bessel functions satisfy orthogonality in form ^U] 

7 

/ rJi{a ln >r)Ji{ai n r)dr = — [jf{ai n r ) - Ji-i(ai n r Q ) Ji+i(ai n raj] 5 n > n = N„iS n > 
o 

(23) 

taking into account boundary condition for optical waveguide. 

Let us now exploit orthogonal relation, we can show that the equation for z 
coordinate are |Hj 

r 2tt 

(D z + al)Al = ^^ [ [rJ l (a ln r)e-^^J2^i m E k E l E m d i pdr, (24) 
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where D z is defined by 



d 2 d 2 

= ^e Q e— 2 - — 2 . (25) 



We choose relation between A\ n and T\ n in form 

d t F? n = ±id z Al, (26) 
this yield to the two orthogonal polarization which can be write as 

E±(r t <p,z,t) = T^J2~ \ d ' A t Jl±i(*lnr)] e # e"+«„ (27a) 
E±(r, <p,z,t) = [d z Af n J l±1 (a ln r)) ^^^ kz + +c.c. (27b) 

l,n 

In this case equations 1127f t have simple form in cartesian co-ordinate system 
Et(x,y,z,t) = t\J2— ^A ± „^±i(«^)e 4(Z±1)v ^ t -^+c.c.,(28a) 

E±(x,y,z,t) = ^J2— d ^i J i±Mnr)e^ l±1) ' p e iut - ikz +c.c., (28b) 



ah 

and we can use it to calculate P z . 

We take into computation electric held in form 

E x = X -Et + \e~, (29a) 
E v = \e+ + \e~. (29b) 

For more generality considering about birefringent axis see 

Let us construct the only transversal mode with fixed a and (3 which mean 
that we chose the simplest form by fixing I — and n = 1 (we cut series to one 
term, this allow us to simplify calculation of the P z ). We introduce a slowly 
varying amplitude of the wave envelope |11| in form 

aX ± ( T ,Oe~ ikz , (30) 

where 

£ = <tz, (31a) 

r = (t-k'z)e, (31b) 

where er is nonlinearity parameter and e is dispersion parameter. 

We insert equations l|28l) into (|30|1 and then we put the slowly varying am- 
plitude l|24|) in it and if the relation between parameters is e 2 ~ a and if we 
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use a new coordinate system which moves at group velocity (l3*T|) . than to the 
second order in e we can obtain nonlinear Schrodinger equation in the form (we 
don't get complex conjugate part) 



tdtX+ + ——d TT X+ + idtX- + ——d TT X- e^M-M = 



7 2 7 »2 



rJ (a ir)—P z (X+,X-)d^dr. (32) 



<r**N 01 k J J ■-»^w at a 
o o 

Let us now evaluate a right site of equations l|32|l . If we keep terms up to the 
order of the third power of e and additional save expression with d T X^ term, 
we have 

id ( X+ ~ik'+d T X+ + ^-d TT X+ + P[\X+\ 2 + Q\X-\ 2 ]X+ = 0,(33a) 
e 2 k" 

id^X- -ik' N d T X- + ^—d TT X- +¥[\X-\ 2 + Q\X + \ 2 ]X- = 0,(33b) 
where 



2 T ° 



&N ik J ' V4" uv V± ' J ' 3a 2 



ijo(a ir) + — j- J 2 (a 01 r) Jq (a ir)fc 2 ] dr(34) 



o 



01 / 



ro 

| r QjoWr)) dr, (35) 
o 



where XPM coefficient (in this case self phase modulation (SPM) equal 1) 

_ 1 3tJ fJiO^OXxxxx^ / „ ( 2 

~~ P 8A^oifc 
and depends on amplitude X*. 

4 Numerical calculation 

First wc define normalized frequency as 



V = - Wei -e 2 (36) 
and coefficient which we calculate 

P = XxxxxCrPcoeff (37) 
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Choosing a value for physical parameters 



w = 12.2 *10 14 Hz (A w 1.54/im) (38a) 
ei = 2.25 (ref. index 1.5) (38b) 
e 2 = 1-96 (ref. index 1.4) (38c) 

e = 8.85 *10~ 12 — (38d) 
m 

N 

Ho = 12.56 *1CT 7 -r-z (38e) 
r from 1.2 * 10~ 6 m to 10 * l(T 6 m (38f) 

On picture n we show results for mode I = and n = (known as TEoi) also 
we show results for mode I — ±1 and n = 1 (known as HEn) which can be 
calculate using above procedure (shows here for I — and n = 0). Picture [5] 
show XPM coefficient. 



mode cut-off line Next mode cut-off line 




Figure 1: Numerical results for Pcoeff for parameters l)38|l. Vertical line shows 
mode cut-off. 

If we consider only one mode with one polarization we can write Nonlinear 
Schrodinger Equation (NLS) 

id z U+^d tt U + 1 \U\ 2 X = 0, (39) 

where 7 is nonlinear coefficient. We compare 7 with work [51 111! IT2) where it 
was defined as 

(40) 
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Figure 2: Numerical results for XPM coefficient l|35|) 



where g depend on variation of the electric field in the fiber cross section and 
in most cases takes a value of approximation 1/2 and 712 is defined by jll) 

3 

712 = 4~ Xxxxx > ( 41 ) 
In our case 7 is define by (for I — ±1 and n = 1 with one polarization) 
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167Vnfciic 2 



We don't have dependence on mode cross-section (radius of light beam and 
fiber) but it is included by the boundary condition. On picture |31 we make 
comparison our results to l|4U|l . On graph G is defined as 

„ 7(our numerical results)c 

G = — — (43) 

gLun,2 



5 Conclusion 

In this paper we shows a new approach to derive a formula for CNLS equations. 

The main idea is to take into account multi-mode case. Here we show the 
simplest case for / = 0, ±1 and n = 1 (as single mode) but if we use equation l|24|l 
and take electromagnetic field (|28l) with more modes we could compute multi- 
mode case. The proceeding is the same as for single but it is more intricate 
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1 2 3 4 5 6 7 

V 



Figure 3: Comparison between results 

and derive simple formula is more difficult (because we have additional terms 
corresponding to mode interaction). 

In this paper we don't show formula for V^ N but it can be calculate. 

Additionally we can allow for change birefringent axis (see eq. Ij29[) ) and 
take case with different grup velocity (introduce fc+ 7^ fc_ ) . 

The definition of the nonlinear coefficient is used in works relating to quan- 
tum effects in fibers |13| , and the coefficient is defined likely in the non-quantum 
NLS equation. 
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